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ABSTRACT 

Gamma-ray binaries emit most of their radiated power beyond ^10 MeV. The non- 
thermal emission is thought to arise from the interaction of the relativistic wind of 
a rotation-powered pulsar with the stellar wind of its massive (O or Be) companion 
star. A powerful pulsar creates an extended cavity, filled with relativistic electrons, 
in the radiatively-driven wind of the massive star. As a result, the observed P Cyg 
profiles of UV resonant lines from the stellar wind should be different from those of 
single massive stars. 

We propose to use UV emission lines to detect and constrain the colliding wind 
region in gamma-ray binaries. We compute the expected orbital variability of P Cyg 
profiles depending upon the interaction geometry (set by the ratio of momentum fluxes 
from the winds) and the line-of-sight to the system. We predict little or no variability 
for the case of LS 5039 and PSR B1259-63, in agreement with currently available HST 
observations of LS 5039. However, variability between superior and inferior conjunction 
is expected in the case of LS 1+61 303. 
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1 INTRODUCTION 

There are now five gamma-ray binaries detected in high 
(HE, > 100 MeV) and very high energy (VHE, > 100 GeV) 
gamma rays: PSR B1259-63, LS 5039, LS H-61 303, HESS 
J0632+057 and IFGL 1018.6-5856. AH contain a massive 
star and a compact object, and emit most of their power at 
energies above 10 M eV (for recent reviews see e.g. iParedesI 
l201ll : lHill et al.ll2010l ). The observed orbitally modulated HE 
and VHE gamma-ray emission indicate that particles are ac- 
celerated to muIti-TeV energies within or close to the binary. 
The available observations in radio. X-rays and gamma-rays 
seem to favor a scenario in which particles arc accelerate d in 
the w ind of a young pulsar (Maraschi & Treves 1981: Dubua 



200E ) rather than in the je t of a microquasar omero et al . 
2OO5I : iParedes et all [2OO6I : iDermer fc Bottcheij I2OO6I ). The 



pul sar model is known to be operating i n case of PSR B1259- 
63 l|Tavani et al.lll994lKirk et al.lll999l '). The latest searches 
for radio pulsat ions in LS 5039 and L S I -1-61 303 yielded no 
positive result l|McSwain et al.ll201ll ). thus the presence of a 
wind collision region and pulsar remains to be proven in all 
remaining gamma-ray binaries. 
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A common denominator to all known gamma-ray bi- 
naries is a massive companion star (O or Be type). The 
massive star loses mass via a strong, hot, radiatively-driven 
wind at a rate that ca n be higher than 10"^ M© yr~^ 
jHowarth fc Priniall989l 'l. UV resonant spectral lines formed 
in the wind are indicators of its velocity and mass-loss rate. 
If the wind density is large enough, the lines have a P Cygni 
profile consisting of a violet-shifted absorption component 
and red-shifted emission component. The P Cyg profiles 
form via absorption and scattering of photospheric contin- 
uum emission in the expanding, spherically symmetric stel- 
lar wind (see Sect. 12. l|) . 

The P Cyg line profiles from massive stars in a binary 
system with another massive star (e.g., 0-f O or WR+0) dif- 
fer from those of a single star and are variable with orbital 
phase. The difference is caused by dynamical interactions 
between the two stellar winds and radiative feedback (e.g., 
wind ionization, radiative braking) between stellar winds 
and radiation field of the companion star. The variability 
of P Cyg pr ofiles in massiv e star binaries was modeled nu- 
merically bv lStevensl(|l993l ). The results showed that P Cyg 
profiles have high diagnostic value for learning on the exis- 
tence and the parameters of the colliding winds region. 

UV line variability is also expected in high-mass X- 
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Figure 1. Section through a generic gamma-ray binary, cut along the binary axis and perpendicular to the orbital plane. Observer is to 
the right at infinity and sees the binary with an inclination i ~ 60°. Solid curves represent surfaces of constant wind velocity component 
towards the observer. The velocities in units of t^z/fcxD are marked below each curve. Negative velocities are in blue, positive in red and 
green is zero velocity. Light gray shaded area marks the extent of the absorption zone. Dark grey area is the occulted wind region, which 
does not contribute to the final profile. Double solid line corresponds to the CD with r] < I. The area to the right of the CD, is the pulsar 
wind cavity, filled with relativistic electrons which do not take part in resonant scattering of UV photons. The yellow ellipse is a cross 
section of one of the emission surfaces Sv . 



ray binaries where a X-ray source ionizes the surrounding 
wind regions creating an extended Stromgren zone. The 
fully ionized plasma does not have any spectral transition, 
therefore the zone does not take part in resonant scatter- 
ing, and the observed P Cyg line profiles should be dif- 
ferent th an in case of single stars. This effect was pre- 
dic ted bv IHatchett fc McCravl (IiOTTI ) and studied in detail 
bv lvan Loon et al.l (|2o'oi[ k 

We propose that line profile variability could be ob- 
served in the UV spectra of gamma-ray binaries within the 
framework of the pulsar wind scenario, where the relativistic 
wind of the pulsar creates a cavity in the companion stel- 
lar wind (see Sect. 12.2")) . The cavity is filled with relativistic 
particles which do not take part in resonant scattering of 
UV photons. Only the wind outside of the cavity is respon- 
sible for P Cyg profile. The cavity is described by its size 
(which increases with the parameter rj - see Sect. 12.21 - and 
thus opening angle), and by location with respect to the 
line of sight. With an increasing size of the cavity, line pro- 
files are expected to differ more from that of a single star. 
Concurrently, because the jets constitute narrow coUimated 
outflows from the vicinity of the compact object, we do not 
expect the creation of the extensive cavity or signiflcant wind 
alterations in the microquasar scenario. The presence or ab- 
sence of line profile variability puts constraints on a size and 
location of the cavity. If changes in the line profile along the 
orbit are detected, a microquasar scenario can be rejected 
with high certainty. 

In this article we study if and how the observed line 
profiles can vary in gamma-ray binaries in the framework 
of pulsar wind scenario. We constructed a semi analytical 



model composed of a massive star wind with a cavity. We 
describe the model in Sect. [2] In Sect. [3] we show how the 
line profiles change with the shape and location of the cavity. 
In Sect. 13.31 we apply the results to gamma-ray binaries LS 
5039, PSR B1259-63 and LS 1+61 303. In Sect.gjwe discuss 
our results and present the conclusions. 



2 MODEL 

2.1 Formation of P Cygni profiles 

In most general terms, the P Cyg line profiles form via line 
scattering of photospheric photons in the hot expanding stel- 
lar wind of a massive star. The absorption line is formed 
when the photons emitted from a photosphere are scattered 
away from the line of sight by the stellar wind located in 
front of stellar disk and moving towards the observer (light 
grey area of Fig.[T]). The absorption component extends be- 
tween Doppler velocities —Voo and 0. The emission line is 
created by the wind which surrounds the stellar disk (as seen 
by the observer) and which scatters the photospheric radia- 
tion into the line of sight, creating an emission line. If there 
was no wind, this radiation would never reach the observer. 
The emission line extends between Doppler velocities —Voo 
and Voo- The radiation scattered behind the star (dark grey 
region of Fig. [T|, does not reach the observer. A superposi- 
tion of the absorption and emission component creates the 
characteristic observed P Cyg profile (Lamcrs & CassineU| 
Il999l ). 

In more detail, the photons of frequency v are scattered 
into or of the line of sight by wind elements moving with 
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respect to the observer with a velocity Vz = v{r)cos6 = 
c(i/ — i'o)/i'o where uo is a rest frame frequency of a given 
resonant atomic transition. All the wind elements moving 
with the same velocity Vz, are distributed along a surface 
Sv The shape and emission of each Si, are symmetrical with 
respect to the line of sight. A cross section across several of 
such constant velocity surfaces, corresponding to different 
frequencies, is shown in Fig. [T] as solid, blue, red and green 
lines. 

A cavity created by a pulsar wind may crop some of the 
constant velocity surfaces, breaking the axial symmetry and 
changing the flux received by the observer from the surfaces. 
The level of modification depends on the shape and location 
of the pulsar wind cavity with respect to the observer. 

2.2 The shape and location of a cavity 

The cavity is bounded by the colliding wind region, com- 
posed of two termination shocks separated by a tangential 
contact discontinuity (CD, double solid line in Fig. [!}. In 
the following we ignore the presence of the shocked winds 
and identify the shape of the pulsar wind cavity with the 
shape of CD. The shape of the cavity can be approximated 
by a conical surface symmetric with respect to the binary 
axis. Its vertex is located at the stagnation point, where the 
ram pressures of the two winds balance: 



E 



AttRI 47rc(s - i?s)2 



(1) 



where Rs is the distance of the stagnation point from the star 
center, and Pp are (respectively) the stellar and pulsar 
wind ram pressures at the stagnation point, s is the binary 
orbital separation at a given orbital phase, M is the stellar 
wind mass loss rate, E is the pulsar spin-down power and 
= v{Rs) is the stellar wind velocity at the radial distance 
Rb from the st ar's centre. For ma ssive star it is given by a 
/3- velocity law (jCastor et al.lll975l ) 



j(r) = 



R*ro 



(2) 



where tq — 1 — {vo/vacY^'^ , vo is the initial wind velocity, 
Vaa is the wind terminal velocity, _R* is the stellar radius and 
/3« 1. 

The size of the cavity can be parametrized by the 
asymptotic opening angle of the cavity t/i, which in the con- 
ical approximation equals to half the vertex angle. We use 
the phenomen ological relation d e rived from hydrodynamical 
simulations bv lBogovalov et al.l (|2008l ). 

7/. = 28.6(4 - r;^/'^)??^/^ for 77^!, 

7/; = 171.5 - 28.6(4 - J7"^''^)ry"^''^ for 77^1, (3) 
where is measured in degrees and where 



V = 



E 



(4) 



is the ratio of pulsar to stellar wind ram pressures at Rs- 
When ri < 1 then tp < 90°, the stellar wind dominates over 
the pulsar wind and the wind-free cavity is created around 
the pulsar as in Fig. [1] When rj > 1 then > 90° , the pul- 
sar wind dominates over the stellar wind and the conical CD 



encloses the star, efficiently blocking the wind from expand- 
ing into space around the binary. In this case the volume of 
the cavity can be larger than the volume of the expanding 
stellar wind. For low enough rj, the stellar wind may also be 
prevented fr om freely expanding ba ckwards by a reconfine- 
ment shock l|Bogovalov et al.1l2008h . Numerical simulations 
are then required to model the interaction region properly 
- this is not taken into account here. Finally, no stable bal- 
ance is possible when rj becomes very high: the pulsar wind 
then rams into the stellar surface (Harding & Gaissor 199(i), 
quenching the stellar wind over at least part of one hemi- 
sphere. 

For a given orbital phase, the location of the pulsar wind 
cavity can be parametrized using only Rs and the angle a 
between the line of sight and the binary axis 



a{i, t, Lj) — cos ^ [sin i cos {t — oj)] . 



(5) 



where t is the pulsar true anomaly (pulsar angular distance 
from periastron) and u is the angle of the line of nodes of the 
orbit. Small values of a correspond to large inclination an- 
gles and orbital phases close to inferior conjunction (INFC, 
when the pulsar is in between the massive massive star and 
the observer). The maximum value is reached at superior 
conjunction (SUPC, when the pulsar is behind the star as 
seen by the observer) with a — 180° when i — 90°. 

2.3 Line profile calculations 

We followed ILucvI (119711 ) (hereafter L71) to compute the 
line profile from the stellar wind of a massive star with a 
radius 7?* and a stellar wind terminal velocity Voo ■ We model 
only single line profiles and do not consider line doublets 
where multi-line scattering takes place. We assumed that 
the resonant transition is possible in the entire volume of 
the wind, as opposed to the situation where line transition 
is only possible in a certain wind layer due to the ionization 
gradient. We assume that the line opacity changes only due 
to changes in the gas density and do not take ionization 
effects into account. Here, we model line profile variability 
only for saturated profiles created in optically thick wind, 
but the qualitative results are also applicable to unsaturated 
lines from optically thin winds. 

To account for the presence of the pulsar wind cavity in 
the stellar wind, we modified the equation 43 of L71 which 
in its original form returns a normalized flux at frequency 
v emitted by a corresponding constant velocity surface Su- 
The modified equation 43 has the form: 

2 r°° 



r dvz , 
X —r^-^rdr 

v(r) oz 



(6) 



where we assumed that the continuum flux equals 1, /i = 
cos 6, H is the Heaviside function, $ is given by equation 
44 in L71, whereas the last term by equation 42 in L71. k 
is an absorption coefflcient set to an arbitrary value which 
gives saturated line proflles. The integral in Eq.[S] includes a 
function which is a fraction of a circular cross section of 
5*1, (at r perpendicular to the line of sight), that is outside 
of the pulsar wind cavity. An example cross section of Sv for 
which = — 0.7foo is shown as yellow ellipse in Fig.[T] The 



© 2009 RAS, MNRAS 000, [l]-?? 



4 A. Szostek, G. Dubus and M. V. McSwain 



Rs [R*] 

2 4 6 8 10 




4 8 12 16 20 

s [R*] 



Figure 2. Contour plot of the normalized profile area Aq as a 
function of the binary separation s and the angle a (between the 
cavity axis and line of sight). Here, = 1 is fixed, which implies 
= 90° and Hs = 0.5s are also fixed. Line profile increasingly 
resembles that of a single star (ylo = l) with increasing s and a. 

function fi^ (r) depends on the shape of the cavity via r], and 
its location with respect to the observer via and a. fu (r) 
is calculated numerically for each binary configuration and 
frequency. 

In order to quantify the difference between P Cyg profile 
P(v) of the same resonant line vq but for winds with different 
cavities, we introduce the profile area A, defined as: 

\P{y) - l\dv (7) 

which is simply a sum of fiux within the emission line above 
continuum and flux removed from continuum by absorption. 
In the following we will use Aq to denote a profile area nor- 
malized by a profile area of a single star. 

A spectral resolution of COS instrument onboard of 
Hubble Space Telescope in the 900 - 1450 A wavelength 
range and using bright object aperture, is i? = 3200 — 4200 
l|Osterman et al.|[201ll '). This corresponds to 40-90 km s~^ 
(2-6% of Voo depending on spectral type). Assuming that 
the line is fully saturated, we estimate that for this resolu- 
tion, a change in profile can be detected if a ratio of profiles' 
areas changes at least by 0.05-0.1. 



3 RESULTS 

3.1 Parameter survey 

We chose to study the variations in line profiles by changing 
T), a and s (or Rs). A. set of three parameters is sufficient to 
describe all possible configurations. We successively freeze 
one of the parameters and show contour plots of Ao as a 
function of the other two. Note that the physical length scale 
is set by 

Figure [2] shows Aq as a function of s and a when r) = 1, 




-1.5 -1. -0.5 0.5 1. 

V [Vco] 



Figure 3. P Cyg profiles for different binary configurations. Top: 
77 = 1, Q = are fixed and the binary separation s is changed. 
Profiles shown s = oo (single star, black, Aq = 1), s = (red, 
Ao = 0.6), s = 6/?* (magenta, 0.7), s = 8R* (blue, 0.8) and 
s = 13-R* (green, 0.9). Bottom: Rs = 3-R*, a = are fixed and 77 
is changed (equivalent to varying cavity opening angle 4>)- Profiles 
correspond to r; = 0.02, ip = 30° (green, Aq = 0.72); r) = 0.23, 
60° (blue, 0.68); r; = 1, 90° (violet, 0.65); -q = 4.5, 120° (pink, 
0.65); 7) = 50, 150° (red, 0.5). 



s [R*] 




Figure 4. Contour plot of the normalized profile area Aq as a 
function of r) and a with Ra = 2/?* fixed. The corresponding 
binary separation s is plotted at the top of the figure. The red 
dashed line shows a = tp. 
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Figure 5. Contour plot of the normalized profile area Aq as a 
function of binary separation s and r] for a = 0° fixed. Constant 
values of Rg in this parameter space are shown as dashed lines. 
The white region at the top of the plot is where the wind in- 
teraction region collides with the stellar surface. The expected 
parameter range for LS 5039 and LS 1+61 303 is also shown. 



which implies — 90° and Rs = 0.5s. The line profile from 
a gamma-ray binary differs more from that of a single star 
when the cavity is closer to the star (smaller Rs in units of R^ 
or, equivalently, smaller s, implying shorter orbital periods) 
and when the axis of the cavity is close to the line of sight 
(low a). Example profiles for q = are shown in Fig. [3^. 
When the cavity moves in closer to the star, the stellar wind 
is prevented from reaching high velocities in the absorption 
zone (Fig. 1). Photons that would have been scattered away 
from the line of sight can then reach the observer unaltered, 
which explains the narrower absorption trough. For rj — 1, 
the red emission component of the line profile barely changes 
since it arises from regions unaffected by the cavity. 

Figure |4] shows Aq as a function of rj and a when 
Rs = 2_R^. Note that keeping Ra fixed implies that the pul- 
sar power increases with t] (hence s), up to _E ~ 10'^* erg 
s~^ when = 10 and assuming M — 10~^ Mq yr^^. The 
line profile is increasingly modified when 77 becomes high (= 
wider cavity opening angle tp). The line profile is most sen- 
sitive to the presence of the cavity when a < i.e. when 
the line-of-sight passes through the cavity for some orbital 
phases. For a > tp, the profile is only weakly sensitive to 
changes of a or -q. When a < 4', changes in the blue absorp- 
tion line (created by the wind located along the line-of-sight) 
are responsible for the strong variations in Aq. Example pro- 
files are shown in Fig. [Sb with Rs — 3i?* and a — 0° . The 
red emission part of the line is modified when 1/) (hence rj) is 
very high, disrupting the wind well beyond the absorption 
zone. 

Finally, Fig. [5] shows Aq as a function of rj and s when 
Q is held fixed at zero. The general trends are in agreement 
with those shown in Fig. [2J11 Aq decreases with increasing 
r] and decreasing s. The dependence on the separation s is 
relatively weak. Instead, the Aq contours follow closely the 
contours of constant Rs, as expected based on the Fig. (2] 



Fig. [5] shows that deviations from a single star profile can 
be observed when Rs < 6i?*, assuming the most favorable 
binary configuration {a — 0°). The white area at the top of 
the plot covers the parameter region where the wind collision 
region intercepts the stellar surface (see §2.2). 

These r esults are in ver y good agreement with the cal- 
culations of IStevenj ()l993l l for binaries composed of two 
massive stars, in particular with t heir Mod e l lb. The pro- 
file variability shown in Fig. 2b of lStevenj l|l993 h. can be 
directly compared with our Fig. [3^. Our calculations apply 
to any inclination, orbital phase and shape of the cavity. 



3.2 Survey limitations 

The li ne profile stron gly depends on the cavity opening an- 
gle -i/". iGavlevI l|2009l ) derived a relation between r/ and the 
asymptotic value of tj), that is different from that used here, 
based on other assumptions about the fiow in the shocked re- 
gion. The opening angles under these assumptions are larger 
from the ones given by Eq. [S] The difference is about 10° for 
0.1 < 77 < 1 and increases to ~ 15° for smaller r; — 0.005. 
Based on Fig. U we estimate that in the most favorable case, 
when a = 0°, the 15° increase of ■(/) leads to a decrease in 
Aq of ~ 0.1, hence of minors consequences to our results. 
Improved accuracy in predicting line profiles is more likely 
to require advanced stellar wind models and numerical sim- 
ulations of the pulsar wind interaction with stellar wind. 

In our calculations we assumed a spherical symmetry 
of the pulsar wind momentum distribution. The simulations 
of pulsar winds from high velocity neutron stars interacting 
with the ISM provide insight on the possible impact of a 
latitudinal variation of pulsar wind on the shape of the in- 
ter action surface. The r esults of the simulations performed 
bv lVieeliusetlD (|20o3) show, that while some changes in 
the structure of the interaction region close to the pulsar are 
possible, the overall geometry remains mostly the same and 
is determined principally by E/c. The present study is far 
from the level of detail where it would be necessary to take 
these effects into account. 

Three out of five known gamma-ray binaries, contain 
a massive Be companion star. The class of stars is known 
for their highly anisotropic winds composed of a fast, tenu- 
ous polar wind and a dense Keplerian equatorial disk. The 
interaction of the pulsar wind with the equatorial disk can 
cause strong variations in momentu m ratio r] and distur - 
bance of the global interaction surface (jOkazaki et akllioilh . 
The shape of the pulsar wind cavity during disk passage and 
its impact on the UV lines are uncertain and can not be ad- 
dressed by our model. 



3.3 Application to known systems 

We apply our model to the gamma-ray binaries with the 
best known parameters: PSR B1259-63, LS 5039 and LS 
1+61 303. We compare line profiles at conjunctions (i.e. at 
minimum and maximum value of a for any given orbit) in 
order to test whether a pulsar wind cavity could be detected. 
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3.3.1 PSR B 1259-63 

This is the only gamma-ray binary known to contain a pulsar 
(spin-down power E — 9, x erg s~^). The companion 
is a Be star, which has a fast, tenuous polar wind and a 
slow, dense equatorial disk wind. It is the polar wind that 
is responsible for the P Cyg profile. 

The polar wind mass-loss rate i n PSR B 1259-63 i s about 
yr~^, fol lowing I Vink et al. 



Mr. 



the stellar parameters from iNegueruela et al. 



corresponding value of r? is ~ 10~ . For this value of ri and 
for i < 25° (Ncgucrucl a et al.ll201ll '). a > ip and we do not 
expect to see changes in line profile (see Fig. [4]and §3.1). 



3.3.2 LS 5039 

UV observations of LS 5039 at two orbital phases, 0.41 (close 
to INFC) and 0.63 (close to apastro n), did not reveal an y 
changes in the P Cyg line profiles (|McSwain et al.l l2004h . 
We set constraints on the binary parameters from the lack 
of variability. 

The value of rj depends on the balance between the 
stellar and (putative) pulsar win d, but the parame t ers o f 
the latter are unknown in LS 5039. ISzostek &: Dubii3 (|201ll ) 
estimated a majcimum value of »7max = 0.6 at periastron 
that, based on Eq. [1] and [3l corresponds to ?/j ~ 75° and 
Rb = 1.2_R*. For higher values of i], the pulsar wind im- 
pinges directly on the stellar surface, disrupting the stellar 
wind over much the pulsar-facing side (the pulsar is only 1-3 
stellar radii away from the surface) - which is excluded by 
the observations. Wit h M = 10"'^ Mp yr~ \ = 2.4 x 10* 
cm s"\ Ri, = 9.3Rq l|Casares et al.ll2005l ) and vo = 2 x 10^ 
cm s~^, ??max is reached for a pulsar with i5 ~ 5 x 10^® 
erg s~^, which is higher than in PS R B1259-63 but not 
unkn o wn among oth e r pul sars (e.g., iKargaltsev fc Pavlov! 
I2OIOI ). IZabalza et al] (|201ll ) set a comparable upper limit 
from the absence of strong thermal X-ray emission from 
the shocked stellar wind. LS 5039 has a mildly eccentric 
orbit, meaning that 77 and i/) decrease and Rs increases at 
phases away from periastron. Hence, line profiles calculated 
with 77 = ?7max at periastron give an upper limit on the 
orbital variability that can be expected. In this case, the 
cavity parameters at conjunctions are (SUPC) -q = 0.38, 
tp = 70°, iis = 1.4i?* and (INFC) tj = 0.17, ip = 55° 
and i?s = 2.8-R*. The inclination of the binary is unknown: 
i < 60° if the X-ray source is point-like (,Casares et al.i,2005, ) 
but i can be as high as 90° if the X-ray source is extended 
l|Szostek fc Dubuil201ll '). However, the observed gamma-ray 
luminosity requires a pulsar spin-down power at least similar 
to PSR B1259-63, in which case rj ^ rj^in ~ 0.04. 

The range of 77 discussed above is marked in Fig. [5l 
showing the expected change in line profile is modest even 
for the most favorable case a = 0°. The difference in P 
Cyg profiles (Fig. [6] between SUPC and INFC for 77 = 77max 
should be detectable down to i ~ 50°. From this result and 
the fact that line profiles of LS 5039 are constant as observed 
by HST, we conclude that 77 must be much less than 77niax 
or that i < 50°. There is a degeneracy between the two 
parameters but the higher the value of i, the lower 77 and 
E need to be to avoid line variability. At high inclination, 
even small cavities with 77 ~ 10"'^ cause an observable {Aq ^ 
0.9) modulation of lines profiles along the orbit. If 77 = 77min 
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Figure 6. Expected P Cyg profiles for LS 5039 and LS 1+61 303 
at two orbital phases (SUPC, red and INFC, blue) assuming r} = 
Tymax (see text for details) and inclination angles of 50° (solid), 
60° (dashed), 90° (dotted). 



then variability in the profiles is detectable unless i ^ 70°. 
Therefore, we conclude that the absence of variability in the 
line profiles observed with HST requires low values of 77 if 
the inclination is high, with an upper limit i ^ 70°, with 
higher values of 77 possible if the inclination is smaller. 



3.3.3 LS 1+61 303 



LS I-I- 61 303 also harbors a massive Be star. I Romero et al.l 
(|2007l ) suggested that the pulsar wind would dominate over 
the polar wind (77 > 1) over much of the orbit, creating 
an extended cavity in the stellar wind. Because the polar 
wind is less dense (lower M) than the wind in LS 5039, 
the line profiles are expected to be weaker but observable. 
The lUE spectra of LS I-I-61 303 (generally low dispersion 
and poor S/N) obtained between 1978-1983 indicate that 
several emission li nes formed in t he wind show variability of 
uncertain nature (|Howarthlll9"8^ ). 

Taking the stellar parameters (7?* = lOM©, i;oo = 10* 
cm s~^, M — 10~*M(7) yrZli orbital solution de- 

rived bv lAragona et al] (|2009l ). we estimate 77max ~ 4.5 at 
periastron, for which Rs = 1.3-R* and = 120°. The spin- 
down power implied is i5 ~ 4 x 10^^ erg s~^. The cavity 
parameters change from (SUPC) 77 = 1.5, Rb = 4.4/?* and 
= 100° to (INFC) ?7 = 2.5, Rb = 1.87?* and ip = 110°. The 
expected LS I-I-61 303 profiles at conjunctions are plotted in 
Fig-E] In Fig. Owe also indicate a range of 77 for values of E 
between 8 x 10^^ and 4 x 10^^ erg s~^ at a = 0. We expect 
to see a significant variability of line profiles along the orbit 
in LS I-I-61 303 because 77 is always relatively high in this 
system for plausible values of E. If the pulsar interacts with 
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the Be disc instead of the polar wind at some orbital phases 
then r; is expected to be much lower than 77max due to the 
increased density, resulting in an increase of Aq. This may 
result in some of the unexplained UV variability although 
this is difficult to assess with the currently available data. 



4 CONCLUSIONS 

We modeled the line profile orbital variability due to the 
presence of a pulsar wind-generated cavity in a gamma-ray 
binary. The line profiles increasingly differ from those of a 
single star when the cavity is largest and located closer to 
the line of sight, that is for smaller Rs, larger rj and smaller 
a. The strong dependence on a means that significant vari- 
ations are expected whenever the cavity intersects the wind 
region responsible for the blueshifted absorption component. 
The emission component of the P Cyg profile remains steady 
except for very large values of (and -0). 

We applied our model to the binaries LS 5039, PSR 
B1259-63 and LS 1+61 303. We predict there is no observable 
line profile variability in PSR B1259-63 due to large sepa- 
ration along the orbit and small inclination angle. HST UV 
spectroscopy of LS 5039 did not show a change in line profiles 
between phases close to INFC and apastron (McSw ain et al.l 
|2004| ). This implies that the inclination is lower than 50° for 
the estimated maximum 77max = 0.6 and no larger than 70° 
if r] — 0.04. For lower values of ri the spin-down power of the 
pulsar is insufficient to explain the 7-ray luminosity. 

Our model calculations predict a significant line pro- 
file modulation in LS I-I-61 303 if > 1 as argued by 
iRomero et al](l2007l l. The modulation should be observable 
even if the binary inclination angle is low. The shape of the 
cavity may be drastically changed during disk wind passage, 
which should be visible as the line profile returning to the 
single star profile. Monitoring of line profile variability with 
orbital phase in LS I-I-61 303 could be used not only to learn 
about the existence and parameters of the pulsar wind cav- 
ity, but also to search for signatures of interactions with the 
equatorial wind of the Be companion. In addition, if changes 
in the line profile along the orbit are detected, a microquasar 
scenario can be rejected with high certainty. 
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